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ABSTRACT 


Acoustic  ‘windows’  are  used  in  SONAR  applications  to  absorb  structural  loads 
associated  with  the  platform  operation  while  allowing  the  passage  of  acoustic  signals. 
The  performance  metric  commonly  used  to  gauge  the  acoustic  window  quality  is 
insertion  loss.  This  thesis  provides  a  derivation  of  insertion  loss  for  multi-layered 
materials  as  a  function  of  frequency  and  angle  of  incidence.  Derivations  are  modified  to 
include  attenuation  of  the  signals  in  the  material  and  the  result  has  been  written  into  a 
MATLAB  model.  Measurements  on  single  layer  plastic,  polyurethane  and  steel  panels 
show  good  experimental  agreement  with  the  theoretical  model.  The  model  is  then  used  to 
predict  insertion  loss  of  multiple  layers  as  a  tool  for  improving  window  bending  rigidity. 
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I.  INTRODUCTION 


A.  ACOUSTIC  WINDOWS 

Acoustic  ‘windows’  and  ‘domes’  are  used  in  Naval  and  commercial  applications 
to  protect  associated  underwater  acoustic  transducers  and  hydrophones  from  afloat  debris 
and  hydrodynamic  loads,  as  well  as  to  streamline  the  structure  exposed  to  flow  loads 
thereby  reducing  flow  noise.  Windows  are  mounted  in  the  frame  of  the  hull,  typically 
positioned  a  few  inches  in  front  of  the  transducer/hydrophone  which  is  in  turn  mounted  in 
a  flooded  seachest.  Tradeoffs  must  be  made  between  required  mechanical  properties  that 
support  operational  loads,  and  a  need  for  optimal  acoustic  performance.  These  two 
requirements  are  typically  mutually  exclusive. 

1.  Operational  Loads 

Operational  loads  fall  into  two  categories:  structural  and  environmental. 

Structural  loads  are  typically  considered  in  terms  of  normal  pressures  applied  to 
the  surface  of  the  window.  Such  loads  are  caused  by  either  hydrodynamic  forces 
associated  with  the  speed  of  the  platform,  or  ‘wave  slap’  applied  to  stationary  platforms. 
Although  stress  in  the  material  is  a  consideration,  a  larger  concern  is  the  displacement  or 
strain  of  the  window  caused  by  bending  under  these  loads.  Significant  changes  in 
geometry  caused  by  hydrodynamic  loads,  may  degrade  the  streamlining  of  the  window  or 
dome  shape  which  in  turn  could  lead  to  unintended  load  and  flow  noise  as  well  as  cause 
damage  to  the  window’s  associated  transducer.  Displacements  should  therefore  be 
minimized. 

Environmental  loads  are  in  the  form  of  external  operating  conditions  such  as 
temperature,  hydrostatic  pressure  and  long  term  seawater  immersion.  Windows  must 
operate  typically  over  a  temperature  range  from  0-40°C,  with  structural  stability  required 
over  the  larger  range  of  -40°C  to  60°C.  Many  materials  experience  changes  in  critical 
material  properties,  such  as  strength,  density,  damping  and  elastic  moduli,  as  functions  of 
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temperature.  Further  changes  are  known  to  occur  if  water  is  absorbed  by  the  material. 
Additionally,  some  materials  experience  significant  dimensional  (thickness)  changes 
under  hydrostatic  load.  Successful  designs  must  consider  and  address  all  of  these  factors 
in  optimizing  the  acoustic  properties  of  windows  employed. 

2.  Acoustic  Performance 

Acoustic  performance  of  the  window  has  a  direct  impact  on  the  associated 
SONAR  system  performance.  Ideally,  acoustic  signals  scattered  from  a  target  should  pass 
through  the  transducer  window  unchanged.  In  practice,  signal  strength  varies  as  a 
function  of  the  frequency  and  the  incidence  angle  of  the  signal.  Degradation  of  the  fonner 
reduces  the  range  at  which  targets  can  be  identified  by  reducing  the  signal  to  noise  ratio 
of  the  system,  while  uncertainty  in  the  latter  may  cause  ambiguity  of  target  direction. 
These  effects  can  be  demonstrated  by  considering  a  simple  line  array. 

Ziomek  [1]  derives  the  normalized  far  field  directivity  (Dn)  function  for  a  line 
array  using  the  spatial  Fourier  transform  of  the  complex  aperture  function,  A(f,xa), 
defined  as  the  complex  signal  produced  by  the  array  in  response  to  a  plane  wave  at  the 
spatial  location  xa  and  at  the  frequency  f.  For  the  case  of  spatial  elements  with  identical 
response,  A(f,xa)  is  defined  as  a  “rectangular”  (constant)  amplitude  window  function,  and 
the  spatial  Fourier  transform  produces  a  sine  directivity  function,  i.e., 


Dn  =  sine  (  s  sin  cp  ) 


s  =  array  length/wavelength 
(p  =  angle  of  incidence 
measured  from  the  nonnal 
to  the  interface  surface. 


Equation  1 .  Line  array  directivity  function 


When  an  acoustic  window  is  in  the  signal  path,  the  aperture  function  used  in  the 
derivation  of  the  directivity  function  should  be  modified  as  follows  to  account  for  the 

amplitude  changes  as  a  function  of  incidence  angle: 
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A  (f,  xa)  — >  A  (f,  xa,cp)  =  rect  (xa  /L)  F((p) 

L  =  length  of  aperture 

Equation  2.  Modified  aperture  function 

The  resulting  directivity  function  would  then  be  the  sine  function  of  Equation  1 
scaled  by  F((p).  The  scaling  function  turns  out  to  be  equivalent  to  the  square  root  of  the 
intensity  transmission  coefficient  (  T  )  metric,  defined  as  the  ratio  of  acoustic  intensity 
transmitted  to  the  intensity  that  was  incident  [2],  Figure  1  and  Figure  2  demonstrate  the 
effect  of  scaling  the  sine  function  by,  in  this  case,  the  square  root  of  the  intensity 
transmission  coefficient  of  a  steel  plate  with  an  observed  significant  sidelobe  level 
increase. 


Figure  1.  Intensity  transmission  coefficient  for  steel  plate  (1  cm  thick)  at  250kHz. 
Angles  are  measured  from  normal  to  the  interface. 
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Figure  2.  Directivity  function  for  a  line  array  for  s  =10  positioned  behind  a  steel 

plate  (1  cm  thick)  at  250kHz. 


A  widely  used  metric  to  gauge  window  performance  is  ‘insertion  loss,’  defined  as 
the  logarithmic  ratio  of  sound  intensity  at  a  fixed  location  with  and  without  the  window 
in  the  path  of  the  sound  wave  [3]. 


Insertion  Loss  =  20  logio  [PI/  P2]  =  10  logio  [1/  F] 

P 1  =  Pressure  at  fixed  point  without  window  in  signal  path 
P2  =  Pressure  at  fixed  point  with  window  in  signal  path. 

Equation  3.  Typical  insertion  loss  definition 

Typical  insertion  loss  levels  on  naval  acoustic  windows  is  2dB  or  less.  Using  an 
electrical  impedance  analogy,  a  characteristic  impedance,  defined  as  the  product  of  the 
material  density  (p)  and  the  material  sound  speed  (c),  is  commonly  compared  with  the 
characteristic  impedance  of  seawater  as  an  indicator  of  window  material  performance.  A 
second  indicator  is  the  ratio  of  material  characteristic  dimension  to  the  wavelength. 
Ratios  of  characteristic  dimensions  of  the  window  to  wavelength  that  are  much  less  than 
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one,  produce  near  acoustic  transparency.  As  such,  in  many  ‘low’  frequency  applications, 
materials  with  poor  impedance  match  to  water  are  still  capable  of  overall  satisfactory 
performance. 


B.  TYPICAL  MATERIALS 

Materials  in  common  use  for  acoustic  windows  fall  into  the  categories  of  rubber, 
polyurethane,  fiber  reinforced  epoxy  and  engineering  plastics.  The  first  two  groups, 
rubber  and  polyurethane,  typically  have  good  impedance  (pc)  matches  to  water.  Table  1 
lists  some  common  materials  along  and  their  associated  ratio  of  characteristic  impedance 
to  that  of  water. 


Material 

Density  (p) 
g/cc 

Longitudinal  Sound 
Speed  (c  ) 
m/s 

(pc)  material /  (pc)  water 

Rubber  (Buna-n) 

0.90 

1550 

1.0 

Polyurethane 
(PR  1574) 

1.01 

1700 

1.02 

Plastic 

(Noryl  EN-265) 

1.08 

2103 

1.1 

ALUM  (  6055) 

2.7 

6300 

11 

Steel  (A36) 

7.7 

6100 

30 

Table  1.  Characteristic  impedance  match  of  common  structural  materials  to  water 


As  stated  above,  there  are  frequency  regimes  and/or  optimization  methods 
available  that  allow  the  use  of  ‘poorly’  matched  materials  in  cases  where  structural 
requirements  exclude  the  use  of  materials  with  better  impedance  matches.  Table  2 
demonstrates  the  mutually  exclusive  nature  of  acoustic  and  mechanical/structural 
properties.  Note  that  materials  with  good  matches  (pc  ratio  ~1),  per  Table  1,  also  have 
relatively  poor  tensile  strength  and  modulus  values. 
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Material 

Tensile  Strength 
(ksi) 

Young’s  Modulus1 
(ksi) 

Rubber  (Neoprene  WRT) 

1.5 

0.4 

Polyurethane  (PR  1574) 

2 

0.8 

Plastic  (Noryl  EN-265) 

6 

390 

ALUM  ( 6055) 

40 

10000 

Steel  (A36) 

36 

30000 

Table  2.  Mechanical  properties  of  common  structural  materials 


Windows  with  lateral  dimensions  approaching  a  wavelength  have  the  added 
complication  of  exciting  flexural  waves  in  the  window  material.  Bounds  applied  in  this 
thesis,  described  below,  will  lessen  these  effects  to  the  point  that  they  may  be  neglected. 


Thesis  Bound  # 

Condition 

Condition 

1 

panel  thickness 

<  2 in. 

2 

frequency  range 

50-250kHz 

3 

angle  of  incidence 

0-40° 

Table  3.  Thesis  Bounds 


C.  PREVIOUS  WORK 

Research  related  to  the  general  case  of  sound  propagation  through  layered  media 
is  extensive.  A  detailed  list  of  early  work  is  reported  in  [4].  Specific  examples  include 
work  by  Fay  [5]  who  presented  a  set  of  equations  to  predict  transmission  loss  for  a  mono- 
layer  plate,  including  comparisons  of  predicted  and  measured  values  for  a  steel  plate 
demonstrating  reasonable  agreement.  This  simple  modeling  has  been  sufficient  for  many 
applications  since  mono-layer  structures  of  rubber,  polyurethane  or  plastic  have  been 
shown  to  satisfy  acoustic  and  structural  requirements  in  earlier  window  designs. 

More  recently,  multi-layer  designs  have  become  necessary  to  meet  strength  and 
bending  requirements  unmet  by  mono-layer  design  windows.  These  multi-layered 
windows  are  considered  composites  and  range  from  typical  fiber  reinforced  epoxy 
(‘fiberglass’)  to  custom  designs.  Fiberglass  designs  have  been  acoustically  characterized 

1  Quasi-static  modulus  listed. 
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empirically  for  specific  applications  and  tend  to  be  used  for  very  low  (<1) 
thickness/wavelength  ratios  because  acoustic  perfonnance  of  these  designs  is  otherwise 
poor.  The  nature  of  these  low  thickness/wavelength  ratios  requires  lower  frequency 
transducers  that  are  typically  larger,  leading  to  larger  dimensions  of  window  aperture 
requiring  in  turn  increased  strength  as  load  pressures  need  to  be  applied  over  a  greater 
area.  For  higher  thickness/wavelength  ratio  applications,  customized  windows  such  as 
those  described  in  [6],  where  a  modified  epoxy  resin  is  often  used  with  a  low  density 
fiber  (‘Spectra’)  resulting  in  a  window  with  macroscopic  density  and  sound  speed  that  is 
nearly  ideal,  have  been  used.  Although  these  newer  designs  offer  similar  acoustic  with 
improved  structural  perfonnance  relative  to  mono-layer  designs,  the  ability  to 
manufacture  consistent  quality  customized  composites  has  been  poor.  Recently  Goodrich 
Corp.  [7]  has  designed  a  multi-layered  window,  called  RHO-COR®,  that  utilizes  an 
elastomeric  ‘core’  material  bounded  by  layers  of  composite  material  and  has  been  shown 
to  satisfy  acoustic  and  structural  requirements  for  frequencies  up  to  50  kHz.  A  drawback 
is  the  fact  that  it  is  costly  relative  to  mono-layer  designs. 
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II.  THEORY 


A.  SINGLE  LAYER  NORMAL  INCIDENCE 

In  the  case  of  nonnal  incidence  sound  waves  passing  from  water  through  an 
isotropic,  homogeneous  solid,  all  points  along  the  water-solid  interface  move  in  phase  so 
no  shear  (transverse)  wave  is  generated  and  the  problem  is  significantly  simplified. 

As  presented  in  [3],  the  intensity  transmission  coefficient  (  T  )  for  a  water  bound 
homogeneous  material  is  found  to  be: 

r  |p  n, transmitted  /  P  n,incident  |  1  /  [1  +  ,25(  [  (| p2c2/pici  )  -  (pic,/p2c2)  f  sin2  k2Lt  ] 

where:  to  =  2  7if,  c  /  =  speed  of  sound  in  water 

c?  =  speed  of  sound  in  material,  k?  =  co/c2,  L,=  thickness  of  layer 

Equation  4.  Transmission  coefficient  for  normal  incidence  case  for  a  water  bound 

panel. 

The  transmission  coefficient  is  related  to  the  insertion  loss  as  follows 

Insertion  Loss  =  10  logio  [1  /  T  ]  =  10  logio  [1  +  .25(  \(p2c2/  p/C/  )  -  (p/c/ p2c2)]2  sin2  k2LJ 
Equation  5.  Insertion  Loss  for  normal  incidence  case  of  a  water  bound  panel. 

Inspection  of  Equation  5  reveals  the  indicators  mentioned  above.  Namely,  closely 
matched  materials  cause  the  \(p2c2  /  p jci  )  -  (p/C//  p2c2)  ]  term  to  approach  zero,  and 
likewise,  k2L  products  near  n  n/2  for  n= 0,1  ,2, 3... cause  the  sin  term  to  vanish.  Both  have 
the  effect  of  causing  insertion  loss  to  approach  zero  which  corresponds  to  “perfect” 
transmission  of  the  acoustic  signal.  On  the  other  hand,  the  insertion  loss  goes  to  infinity 
as  r  approaches  zero,  with  the  physical  interpretation  in  this  case  that  little  sound  is 
allowed  to  pass  through  the  material.  Figure  3  demonstrates  this  point. 
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Transmission  Coefficient 


Figure  3.  Insertion  loss  v  transmission  loss 

Such  results  are  useful  in  the  idealized  case  of  pure  normal  plane  wave  incidence. 
The  more  common  case  is  wave  incidence  that  is  not  strictly  normal,  and  in  such 
instances,  a  shear  wave  is  generated  in  the  solid  because  points  along  the  water-solid 
interface  are  not  excited  in  phase.  For  the  off-normal  incidence  case,  a  more  detailed 
approach  that  requires  an  understanding  of  stress  tensors  is  necessary. 

B.  MULTI-LAYER  NON-NORMAL  INCIDENCE 

1.  Stress/Strain  Wave  Propagation 

Waves  in  solids  can  be  described  as  propagations  of  stress  and  strain  within  the 
material.  Officer  [8]  provides  a  detailed  derivation  of  the  governing  wave  equation  in 
solids.  Brekhovskikh  [9]  expands  this  treatment  by  considering  wave  propagation 
through  multi-layered  solids.  In  both  references,  fundamental  concepts  of  stress  and 
strain  are  required  as  groundwork  in  the  equation  derivation.  In  the  following,  a  short 
derivation  of  the  governing  equations  is  given. 

Stress  (T)  is  defined  as  the  force  acting  on  a  given  area  (8A)  divided  by  that  area. 
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T  =  f/8A 

Equation  6.  Stress  (T)  definition 

The  stress  state  at  a  point  in  a  continuum  can  be  resolved  into  three  traction 
vectors  (t1,  t2,  t3),  each  acting  on  a  differential  area  whose  normal  is  one  of  the  three  axis. 


Figure  4.  Traction  vectors  acting  on  orthogonal  planes  of  a  differential  element  from 

[10]. 


Each  of  these  three  tractions  can  then  be  further  resolved  into  three  orthogonal 
components.  The  resulting  nine  stress  components  are  referred  to  using  the  following 
indicial  (tensor)  notation:  T;j,  where  i  =  index  of  the  nonnal  plane,  and  j  is  the  coordinate 
component  of  the  applied  stress. 
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Figure  5.  Traction  vector  components  acting  on  a  differential  element  from  [10]. 


Before  the  equations  of  strain  can  be  fonnulated,  a  definition  for  media 
displacement  is  required.  Therefore,  the  displacement  u,  identifies  the  displacement  of  a 
point  in  the  i  direction.  Extensional  strain  can  generally  be  defined  as  the  relative  change 
in  position  per  unit  distance.  In  a  simple  case,  the  change  in  axial  length  (SuO  divided  by 
the  total  length  (Sx).  Shear  strain  is  defined  differently,  as  shown  in  figure  below  and  in 
Equation  7. 


Shear  strain  tensor  is  the  average 
of  two  strains,  i.e., 

s  12  =  s  21  =  ‘<4  (Su2/Sx  +  6ui/Sy) 

Figure  6.  Shear  strain  definition  from  [10]. 
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In  a  similar  fashion  to  stress,  the  strains  are  referred  to  by  indicial  notation:  s 


sn  =  8ui/8x  s  22  =  8u2/8y  s  33  =  8U3/8Z 

s  12  =  s  21  =  Vi  (8U2/8X  +  8ui/8y ) 

S  13  =  £  31  =  Vi  (Sui/8z  +  SU3/8X  ) 

S  32  =  £  23  =  Vi  (Su3/8y  +  SU2/8Z  ) 

Equation  7.  Strain  notation 

The  generalized  fonn  of  Hooke’s  law  is  used  to  relate  each  of  the  nine  stress 
components  as  a  linear  function  of  the  nine  components  of  strain,  each  multiplied  by  a 
‘stiffness’  coefficient,  cpi.  Officer’s  [8]  derivation  makes  use  of  symmetries  in  stress  and 
strain  tensors  identities  and  equilibrium  considerations  to  reduce  the  generalized  Hooke’s 
law  relations  to  six  equations,  displayed  in  Equation  8.  A  widely  used  enumeration 
scheme  to  reduce  the  number  of  subscripts  on  each  term  is  given  below. 


Tl=  Cl  1  8  1  +  C12  82  +  C13S3+  Ci4  84+  C15  85+  Ci6  £6 
T2=  C21  £  1  +  C22  £2  +  C23£3  +  C24  £4+  C25  £5  +  C26  £6 
T3=  C31  8  1  +  C32  £2  +  C33  S3  +  C34  84+  C35  85+  C36  £6 
T4=  C41  8  1  +  C42  £2  +  C43  S3  +  C44  84+  C45  £5  +  C46S6 
Ts=  C51  8  1+  C52  £2+  C53  £3  +  C54  £4  +  C55  £5  +  C56  £6 
Tg  C61  £  1^""  C62  £2  C63  £3  ^64  £4  ^65  £5  -^-£66  £6 

Equation  8.  Reduced  form  of  generalized  Hooke's  law. 

Nyc  [11]  presents  strain  energy  considerations  that  require  cmn  =  cnm.  When 

applied  to  Equation  8,  this  reduces  the  number  of  independent  stiffness  coefficients  to  21. 

13 


Certain  material  classes  allow  these  equations  to  be  reduced  even  further  if  independence 
of  one  or  more  coordinate  directions  in  the  medium  exists.  The  table  below  lists  some 
material  classes  and  the  number  of  independent  stiffness  coefficients. 


Crystal  Type 

Use 

Material  Example 

Triclinic 

21 

Axinite 

Monoclinic 

13 

Lithium  Sulfate 

Trigonal 

6 

Quartz 

Hexagonal 

5 

Piezoelectric  ceramic 

Transversely  isotropic 

5 

Unidirectional  Layer  in  Composite  Material 

Isotropic 

2 

Steel,  Engineering  Polymers,  Elastomers 

Table  4.  The  number  of  independent  stiffness  coefficients  (nsc)  for  various  crystal 

types. 


Officer  [8]  reduces  Equation  8  to  two  independent  elastic  constants  for  isotropic 
materials;  X  and  p,  known  as  Lame’s  constants. 


Ti=  (2p  +  X)  s  i  +  X  82  +  X  83 
T2=X  s  1  +  (2 p  +  X)  82  +  X  S3 
T3=  X  8  1  +  X  82  +(2p  +  X)  S3 

T4=  p  84 

T5=  p  £5 
T6=  ps6 


Equation  9.  Isotropic  stress-strain  equations. 


2.  Non-Normal,  Multi-Layered,  Unattenuated  Case 

For  the  case  of  a  planar,  liquid-solid  interface,  normal  and  tangential  components 
of  stress  must  be  equal  at  that  interface.  However,  since  fluids  support  shear  stress  to  a 
negligible  degree,  the  tangential  components  of  stress  at  the  boundary  are  set  to  zero 
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while  normal  stress  balances  fluid  pressure.  A  third  boundary  condition  is  continuity  of 
the  normal  components  of  displacement  at  the  interface.  Conditions  are  summarized 
below: 


1 .  Continuity  of  normal  stress  (T3 )  across  the  boundary 

2.  Tangential  stress  (T5 )  of  zero  at  the  boundary. 

3.  Continuity  of  normal  displacement  (u3  ) 


Table  5.  List  of  boundary  conditions  for  general  case  of  nonnal  and  non-normal 

incidence. 


For  two  dimensional  problems  that  are  independent  of  the  y  coordinate,  one  need 
only  consider  the  x-z  plane  of  the  medium.  For  any  angle  of  incidence  as  shown  in  Figure 
7,  Brekhovskikh  derives  amplitudes  of  plane  waves  propagating  through  the  layered 
material  by  repeated  applications  of  the  boundary  conditions  of  Table  5  at  every 
interface. 


Figure  7.  Wave  propagation  direction. 
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Brekhovskikh’s  derivation  at  a  single  interface  leads  to  a  set  of  linear  equations 
presented  in  matrix  form  in  Equation  10.  It  is  important  to  note  that  the  incident  wave 
propagates  in  the  negative  z  direction. 


coT572J 


jM  cos(Pz)  -  M  sin(Pz) 


-P  sin(Pz)  jPcos(Pz) 


Q  sin  (Qz) 


jM  cos(Qz) 


-jQ  cos(Qz)  f  ((j)’ +()>”) 
-  M  sin(Qz)  (t))’  -  (J)”) 


-j  K  cos(Pz)  Ksin(Pz)  2|j  MQ  sin(Qz)  -  j2|jMQ  cos(Qz)  (V  ) 


MPsin(Pz)  -jMP  cos  (Pz)  j  ((Q2-M2)/2)  cos(Qz)  ((M2  -Q2)/2)  sin(Qz) 


(y  -v”) 


Or,  in  short  form: 


r11]  =  [An]  [sn] 

6 

where:  P( 

Mn  =  k„  sinyn  =  kn  sin  cpn=  k„+i  sin  cpn+i  ^ 
Pn  =  kn  cos  cpn  ^ 

Qn=  Kn  cos  y „  7 

K  =  (M2k+P2k+2pP2)  7 

cL  =  Compressional  wave  speed  7 

Cg  =  Shear  wave  speed 
kn  =  co  /cL  in  layer  ‘n’ 

Kn  =  00  /cg  in  layer  ‘n’ 

yn  =  Shear  wave  propagation  angle  in  layer  ‘n’ 
opn  =  Compressional  wave  propagation  angle  in 
layer  ‘n’ 


(|)  =  Compressional  velocity 

potential 

(j)’  =  Incident 

()>”=  Reflected 

\|/  =  Shear  velocity  potential 

\|/’  =  Incident 

\|/”=  Reflected 


Equation  10.  Matrix  form  of  particle  velocities  and  stress  component  equations  in  an 

isotropic  layer,  n. 

A  key  concept  used  by  Brekhovskikh  in  this  multilayered  scenario  is  that  the  [Sn] 
do  not  have  a  z  dependence,  and  therefore  are  constant  through  each  layer.  This  means 
that  stress  at  the  top  of  the  layer  can  be  related  to  the  stress  at  the  bottom  of  a  given  layer 
by  elimination  of  this  common  factor. 


16 


The  A  matrix  in  Equation  12  is  a  4x4  matrix  that  is  the  product  of  [E1]. 

Finally,  for  the  special  case  of  isotropic  material  bounded  by  water,  Brekhovskikh  derives 
a  pressure  ratio,  which  can  be  related  to  the  intensity  transmission  coefficient,  which  in 
turn  can  subsequently  be  used  to  determine  insertion  loss  at  both  normal  and  non-normal 
angles  of  incidence. 
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Insertion  Loss  =10  logio  [1  /  F  ] 


r=[-2Pn  co2pi  /  (H1+H2)  ]2 


where  : 

pi  =  Density  in  layer  1  (seawater) 
Hl=Pn[PiY3  -  co  pi  co  Y4] 
Y1  =  [(-I/A41  )  A42  A21+  A22  ] 
Y3  =  [(- I/A41  )  A42  A31+  A32  ] 


pi  =  Density  in  layer  n 

H2  =  -pnco2  [Pi  Y1  -  co  pi  co  Y2] 

Y2  =  [(-I/A41  )  A43  A21  +  A23] 

Y4  =  [(-1/A4i)  A43  A31  +  A33] 


Equation  13.  Insertion  loss  using  intensity  transmission  coefficient  for  multi-layered, 
water  bound,  isotropic  material  at  normal  and  non-normal  angles  of  incidence. 


3.  Non-Normal,  Multi-Layer,  Attenuated  Case 

Equation  5  and  Equation  13  neglect  attenuation  of  the  pressure  wave  as  it 
propagates  through  the  layer.  Brekhovskikh’s  approach  simply  multiplies  potentials  by 
decay  terms,  such  as  cxp(-cz),  where  2  is  the  attenuation  term  and  z  is  the  spatial 
variable.  The  shear  and  compressional  decay  terms,  and  £,L ,  are  defined  by  Capps  [12] 
below  in  tenns  of  the  loss  properties  in  the  material,  characterized  by  the  ‘tan  8’  metric. 


^G=  SQRT  [  pco2  ((  1  +  tan2  8G  f  -  1)  /  2  (f  (  1  +  tan2  8G  )  ] 

=  SQRT  [  pco2  ((  1  +  tan2  8L  )'A  -  1)  /  2L’(  1  +  tan2  8L  )  ] 

L  =  (A,  +  2p  ) 

L’=  Real  part  of  L 
L”=  Loss  part  of  L 

tan  Sl  =  loss  tangent  associated  with  L  =  L”/L’ 
tan  Sq  =  loss  tangent  associated  with  p  =  p’Vp’ 


Equation  14.  Compression  and  shear  decay  term  definition. 


Since  the  incident  wave  is  traveling  in  the  negative  z  direction,  the  resulting  wave 
propagation  takes  the  following  form,  using  the  variables  from  Equation  10: 
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(j)  =  [f  exp  (jPz  -E,l (d-z»  +  (j)”  exp  (- jPz-^Lz)]  exp  (jM-^Lx) 

(j)  =  [(j)’  exp  (jPa  z  -^L  d)  +  (()”  exp  (- jPa  z)]  exp  (jMa  x) 
where:  Pa  =  P  -j^L  ,  Ma  =  M  -j£,L 

V  =  [v)/’  exp  O'Qz  -^o(d-z))  +  v|/”  exp  (- jQz-^Gz)]  exp  (j'M-^gx  ) 

\|/  =  [\|/’  exp  (jQa  z  -^G  d)  +  \)/”  exp  (-  jQa  z)]  exp  (jMa  x  ) 
where:  Qa  =  Q-j£,G 

Equation  15.  Velocity  potentials  equation  for  the  attenuated  wave  case. 

When  attenuation  exists,  a  modified  set  of  equations,  based  on  the  same  approach 
used  to  develop  Equation  10,  can  be  derived.  Equation  16  displays  these  modified 
equations,  based  on  the  modified  potentials  of  Equation  15,  in  matrix  form.  It  can  be  seen 
that  the  potential  sums  shown  in  Equation  10  are  replaced  by  individual  potentials  for  the 
upward  and  downward  traveling  compressional  and  shear  waves,  scaled  by  decay  terms. 
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V; 


Ma  el 


Ma  e2 


-  Qa  e3 


Qa  e4 


-Ka  e2  -2|_i  MQa  e3 


2|aMQa  e4 


-2)jPa  Ma  el  2|jPa  Ma  e2  -|a  (Ma2  -Qa2)  e3  -  (Ma2  -Qa2)  e4 


Or,  in  short  form: 

[T]  =  j  [A2]  [S2] 


where: 

Ka  =  (Ma\+Pa\+2|uPa2) 


el  =  exp  ( jPaz-£,Ld) 
e2  =  exp  (-jPaz) 
e3  =  exp  (  jQaz-cGd) 
e4  =  exp  (-jQaz) 


Equation  16.  Matrix  fonn  of  linear  equations  with  upward  and  downward  potentials 

separated. 


A  recurrence  relationship  similar  to  that  described  in  Equation  12  also  exists  for 


the  attenuated  case. 
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[T  z=bottom  of  layer  n- 1  ]  [A2  z=  bottom  of  layer  n- 1  ]  [S2  ] 

[S2n_1]  =  [A2n" 

z=  bottom  of  layer  n- 1  z=bottom  of  layer  n- 1  ] 

[T  z=top  of  layer  n-1  ]  =  [A2n_1z  =top  of  layer  n- 1  ]  [S2n_1] 

[T  z=top  of  layer  n- 1  ] 

[A2n“ 

z=top  of  layer  n- 1  ]  [A2n-  z=  bottom  of  layer  n- 1  ]  ''[T11-1  z=bottom  of  layer  n- 1  ] 
[T  z=bottom  of  layer  n  ]  =  [T  z=top  of  layer  n- 1  ] 

[T  z=bottom  of  layer  n  ] 

[A2  z=top  of  layer  n- 1  ]  [A2  z=  bottom  of  layer  n- 1  ]  [T  z=bottom  of  layer  n- 1  ] 

[T  z=bottom  of  layer  n  ]  =  [I'1"1]  [Tn  l  z=bottom  of  layer  n- 1  ] 

[Tn"‘]  =  [In“2]  [Tn"2] 

[Tn]  =  [I”'1]  [Sn“2] . [S1]  [T1] 

=>[Tn]=  [A]  [T1] 


Equation  17.  Matrix  fonn  of  recurrence  relationship  between  isotropic  layers  n  and  1 

for  the  attenuated  case. 


The  intensity  transmission  coefficient  can  be  detennined  using  a  similar  approach, 
modified  to  account  for  the  matrix  differences  between  Equation  10  and  Equation  15,  to 
that  used  by  Brekovskikh  in  the  non-attenuating  case  referenced  above.  Again,  the 
insertion  loss  metric  is  related  to  the  transmission  coefficient.  The  final  equation  is  shown 
below. 
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Insertion  Loss  =10  logio  [1  /  Ta  ] 
ra  =  [  -2  Pa,n  pi  Ka,n  /pn  (H3  +  H4)  f 

where  : 

Ka,n  =  Ka  in  layer  n 
Pa,n  =  Pa  in  layer  n 
Pa,i  =  Pa  in  layer  1 

H3  =  -Ka,n  [-Yli  Pa,i  +  piC02Y2i  ]  H4  =  -  Pa,n  [-Y3i  Pa,i  +  piC02Y4i  ] 

Yli  =  [(-  1/Ai4r  )  Ai  42  Ai2i+  Ai22]  Y2i  —  [(-  1/Ai 41  )  Ai43  Ai2i+  Ai23] 

Y3i  =  [(-  1/Ai 41  )  Ai42  Ai3i+  Ai32]  Y4i  =  [(-1/Ai 41  )  Ai43  Ai3i+  A  i 33] 


Equation  18.  Insertion  loss  with  attenuation  using  intensity  transmission  coefficient  for 
multi-layered,  water  bound,  isotropic  material  at  nonnal  and  non-normal  angles  of 

incidence. 

Since  literature  values  for  X  are  uncommon,  it  can  be  calculated  using  Young’s 
modulus  (E)  and  p.  In  turn,  p  is  found  given  E  and  v  according  to  Equation  19  from 
Loeser  [13],  Due  to  the  lack  of  experimental  data,  a  simplification  is  made  concerning 
Poisson’s  ratio  that  assumes  it  is  a  real  number:  v”  =0  so  v*  =  v’. 


Equation  19.  Relationship  between  shear  modulus,  Young’s  modulus  and  Poisson’s 

ratio  for  the  case  of  v”  =  0. 
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Using  Equation  19  and  the  definition  of  X  in  terms  of  E  and  v  shown  in  Loeser,  X* 
can  be  calculated  as  follows: 


X*  =  v'E*  /  (l+v’)(l-2v’) 

2|u*  =  E*  /  (1+v’) 

(X*  +  2p*  )  =  (E*  /  (1+v’))  [v7  (l-2v’)  +  1] 


Equation  20.  Relationship  between  compressional  modulus,  Young’s  modulus  and 

Poisson’s  ratio  for  the  case  of  v”  =  0. 


Since  signal  decay  is  primarily  associated  with  the  materials’  elastic  moduli, 
introduction  of  complex  terms  affects  the  magnitudes  of  propagation  speeds  in  the 
material.  Capps  defines  these  speeds  using  loss  tangents  as  follows: 


cG  =  SQRT  [  2p’(  1  +  tan2  SG  )  /  p(  (  1  +  tan2  8G  )V2  +  1)  ] 
cL=  SQRT  [  2L’(  1  +  tan2  SL  )  /  p(  (  1  +  tan2  SL  fA  +  1)  ] 


Equation  2 1 .  Sound  speed  equations  for  lossy  material. 


Using  Equation  19  and  Equation  20,  Equation  14  and  Equation  21  can  be 
rewritten  in  terms  of  complex  Young’s  and  Poisson’s  moduli  E  and  v. 


cL=  SQRT  ( [2  (E7  (1+v’)  [v'/  (l-2v’)  +  1])/  p]  [TD/(TD1/2  +  l)] ) 
cG=  SQRT  (  [2/(E’/2  (1+v’))  p]  [TD/(TD1/2  +  l)j  ) 

^G=  SQRT(  [  pco2  12  (E’/2  (1+v’))]  [(TD1/2-l)/TD  ]  ) 

^L=  SQRT(  [  pco2  /2(E’/  (1+v’)  [v7  (l-2v’)  +  1]  )]  [(TD1/2-l)/TD  ]  ) 

TD  =(  1  +  tan2  §E  ) 


Equation  22.  Compressional  and  shear  sound  speeds  and  decay  terms  expressed  as 

functions  of  E  and  v. 


23 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


24 


III.  MEASUREMENTS  AND  MODELS 


A.  MEASUREMENTS 

To  verify  the  theoretical  models  discussed  in  previous  sections,  experimental 
measurements  were  made  on  panels  of  several  materials  listed  in  Table  2.  These 
measurements  were  made  at  several  different  open  water  facilities  using  identical 
measurement  protocols. 

Insertion  loss  was  calculated  according  to  Equation  3  using  sound  pressure  level 
(SPL)  measurements.  An  incident  signal  was  produced  by  a  calibrated  reference  standard 
projector,  with  and  without  a  panel  in  the  acoustic  signal  path  between  the  projector  and  a 
calibrated  reference  standard  hydrophone.  The  rotator  shaft  is  used  to  position  the  panel 
such  that  a  plane  wave  from  the  projector  strikes  the  panel  at  varying  angles  of  incidence. 
The  setup  is  shown  below  in  Figure  8. 
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Figure  8.  Acoustic  Measurement  Setup 
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In  practice,  two  issues  must  be  addressed  to  assure  test  setup  and  panel  sizes  are 
sufficient  to  address  geometric  influences  affecting  measurement  quality.  The  first  of 
these  is  the  fact  that  the  wavefront  incident  on  the  panel  is  not  planar,  but  spherical.  The 
issue  can  be  mitigated  by  applying  the  criterion  that  the  maximum  dimension  of  the 
hydrophone  used  to  measure  the  signal  be  such  that  the  incident  signal  at  the  edge  and  the 
center  of  the  hydrophone  be  no  more  than  30°  out  of  phase,  or  A,acoustic/12.  This 
requirement  forces  a  maximum  hydrophone  length  (Li  max  )  defined  by  Equation  23. 


Equation  23.  Minimum  panel  length  for  plane  wave  assumption. 


The  second  geometric  issue  is  an  effect  called  ‘edge  diffraction’  [3].  Using  Figure  8  and 
Figure  9,  R3  =  Ri-  R2  .  The  distance  from  the  edge  of  the  panel  to  the  center  of  the 

hydrophone  is  shown  in 

Figure  9  as  R3’.  The  panel  must  be  large  enough  such  that  the  difference  between 
R3’  and  R3  is  large  enough  to  allow  a  signal  to  be  sampled  before  edge  diffraction  signal 
arrives.  Since  the  reference  hydrophone  requires  2-3  cycles  to  reach  steady  state,  this 
criterion  is  set  at  four  wavelengths  (minimum)  to  allow  one  full  cycle  to  sample. 
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Figure  10  shows  that  hydrophone  length  should  not  be  larger  than  0.09m  (9cm). 
The  FI-52  type  hydrophone  (Length  =  5cm)  is  the  largest  used  in  these  measurements  and 
easily  meets  the  criteria.  The  line  marked  ‘L2’  in  Figure  10  represents  the  minimum 
panel  length  required  as  a  function  of  frequency.  The  largest  panel  size  required  is  ~.29m 
(11  inches)  at  the  low  end  of  the  band. 


LI_ L2 


Frequency  (Hz) 


Figure  10.  Minimum  required  panel  length  L2  and  maximum  hydrophone  length  L|. 
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Both  transducer  and  hydrophone  have  accuracies  within  +/-  O.ldB,  in  ideal 
conditions.  Two  measurements  are  required  to  establish  the  insertion  loss  value. 
Therefore,  total  inaccuracies  in  measured  data  could  be  as  large  as  +/-  0.4dB. 

B.  MODEL 

Equation  5  ,  Equation  13  and  Equation  18  ,  derived  in  the  previous  section  for 
isotropic  layers  bounded  by  water,  represent  three  levels  of  insertion  loss  models  with 
increasing  levels  of  complexity;  Equation  5  for  the  simplest  case  of  normal  incidence/  no 
attenuation,  Equation  13  for  non-normal  incidence  without  attenuation,  and  Equation  18 
for  non-normal  incidence  that  includes  attenuation.  Each,  hereafter  referred  to  as  Model 
1,  2,  and  3,  respectively,  can  be  used  as  the  basis  equation  for  an  insertion  loss  model  that 
require  user  input  of  the  following  variables: 

1.  Number  of  layers 

2.  Angle  of  Incidence  (Equation  13  and  Equation  18) 

And  for  each  layer: 

3.  Young’s  Modulus  (E) 

4.  Poisson’s  ratio  (v  ) 

5.  Young’s  tan  8  (Equation  18) 

6.  Layer  thickness 

Although  it  is  not  readily  obvious,  a  quick  comparison  between  methods  for  the 
case  of  a  single  layer  of  non-attenuating  steel  that  is  one  cm  thick,  with  a  normally 
incident  acoustic  plane  wave,  produces  identical  results.  The  MATLAB  code  for  model  3 
is  included  in  the  Appendix. 

To  further  validate  Model  3,  it  is  also  checked  against  cases  of  single  layered, 
isotropic  materials  at  various  incidence  angles  and  levels  of  attenuation.  Material 
property  inputs  used  to  test  the  model  are  listed  in  Table  6. 

Density  and  Poisson’s  ratios  are  found  in  multiple  sources  for  the  items  listed  in 
the  table  below.  The  dynamic  properties,  Young’s  (Storage)  modulus  and  tan  8,  vary 
significantly  with  temperature  and  frequency.  Little  data  for  the  specific  materials  and 

frequency  ranges  used  in  this  thesis  can  be  found  in  the  literature.  In  the  case  of 
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polyurethane,  Capps  [12]  provides  tan  5  measurements  versus  a  scaled  frequency  across  a 
broad  frequency  band.  The  scaling  coefficient  accomplishes  a  tan  8  shift  at  temperatures 
away  from  the  reference  temperature.  Specific  moduli  and  tan  8  measurements  on  Noryl 
plastic  were  available  [14].  On  the  other  hand,  nylon  6  moduli  and  tan  8  values  were 
estimated  from  available  data  at  frequencies  outside  the  range  of  the  thesis  bounds. 


Material 

Density(gZcc) 

Young’s 

Modulus 

(GPa) 

Tan  8 

Poisson  ’s 
Ratio 

(v) 

Plastic2 

(Noryl  EN-265) 

1.08 

2.62 

0.029 

0.38 

Plastic3 
(Cast  Nylon  6) 

1.12 

4.14 

0.023 

0.38 

Steel  (A36) 

7.7 

207 

~0 

0.29 

Polyurethane 
(PR  1574) 

1.03 

0.200 

0.150 

0.49 

Table  6.  Physical  properties  for  single  layer,  isotropic  materials  used  in  the  test 

cases. 


C.  NORMAL  INCIDENCE  MODELS  AND  MEASUREMENTS 

In  Figure  11  and  Figure  12,  the  model  results  and  the  measured  data  are  plotted 
for  the  two  plastic  material  cases  subject  to  normal  incidence  acoustic  waves.  The 
agreement  between  the  model  and  the  measured  data  is  good.  The  fact  that  the  measured 
data  shows  a  sinusoidal  shape  of  a  similar  period  as  the  model  data  indicates  the  sound 
speed  and  density  estimates  are  close  to  the  actual.  The  attenuation  values  are  also 
validated,  based  on  the  similar  linear  ‘DC  bias’  evident  in  both  the  model  and  measured 
cases. 


-  Noryl  EN-265  plastic  dynamic  properties  determined  from  sample  measurements  at  various 
temperatures  [14]. 

3  Cast  Nylon  6  dynamic  properties  are  estimates.  Modulus  is  estimated  from  sound  speed 
measurements  made  by  Underwater  Sound  Reference  Division  (USRD)  along  with  density  and  Poisson’s 
ratio.  Tan  5  value  estimate  is  made  using  ultrasonic  (2.25MHz)  measurements. 
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Model  1  - Model  3  —"—Measured  Data 


Frequency  (Hz) 


Figure  1 1 .  Insertion  loss  of  1 .905cm  cast  Nylon  6  panel  at  normal  incidence. 


—  Model  1  —  Model  3  Measured  Data 


Frequency  (Hz) 


Figure  12.  Insertion  loss  of  2.54cm  Noryl  EN-265  panel  at  normal  incidence. 


Normal  incidence  steel  plate  measurements  shown  in  Figure  13  indicate  that, 
although  there  is  fluctuation  in  the  measurement,  the  model  to  measurement  match  is 
good.  This  is  evident  when  the  following  is  considered.  When  insertion  loss  is  high,  as  in 
this  case  of  >  17dB,  the  received  signal  has  an  increased  susceptibility  to  measurement 
system  noise,  both  electronic  and  acoustic.  Electronic  noise  is  related  to  measurement 
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system  equipment  and  signal  to  noise  ratio.  Acoustic  ‘noise’  is  the  result  of  small, 
otherwise  insignificant,  reflections  from  various  scatterers  located,  unavoidably,  within 
such  range  as  to  coherently  sum  with  the  incident  signal.  In  this  case  of  the  steel  plate, 
the  periodicity  of  the  ripple  in  Figure  13  suggests  the  measured  signal  contained  a 
reflection  from  a  nearby  (0.05m)  scatterer,  likely  the  hydrophone  fixture.  If  the  steel 
panel  data  points  are  least  squares  fit  to  a  3ld  order  polynomial,  the  agreement  between 
measurements  and  model  greatly  improves,  as  shown  in  Figure  14. 
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Figure  13.  Insertion  loss  of  1.0cm  A36  carbon  steel  plate  at  normal  incidence. 


31 


Frequency  (Hz) 


Figure  14.  Insertion  loss  of  1.0cm  A36  carbon  steel  plate  at  normal  incidence  with 
measured  data  points  fit  with  a  3rd  order  polynomial  curve. 


Finally,  the  case  of  a  polyurethane  panel  is  presented  below  in  Figure  15.  Note 
that  the  insertion  loss  level  is  very  low,  relative  to  the  steel  plate  case.  The  measured  data 
agree  very  well  with  model  1,  which  does  not  include  an  estimate  on  attenuation  in  the 
material.  Model  3,  on  the  other  hand,  does  include  the  effect  of  wave  attenuation  through 
the  material  based  on  an  attenuation  value  given  by  Capps  [12].  Since  model  1  and  model 
3  inputs,  with  the  exception  of  attenuation,  are  identical,  its  clear  that  the  disagreement 
between  model  3  and  the  measured  data  is  based  on  the  attenuation  input.  If  the 
attenuation  value  is  lowered  by  an  order  of  magnitude,  agreement  between  model  3, 
model  1  and  the  measured  data  greatly  improves,  as  demonstrated  by  Figure  16.  As 
mentioned  previously,  there  is  scarce  measured  attenuation  data  in  the  literature  for  the 
specific  frequency  range  and  material  used  in  this  thesis.  As  such,  the  attenuation  input 
given  by  Capps  and  used  for  the  polyurethane  panel  model  could  not  be  verified  by  a 
second  literature  source  and  its  accuracy,  based  on  the  measurements  done  here,  is 
suspect. 
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Figure  15.  Insertion  loss  of  0.635cm  thick  PR-1574  panel  at  normal  incidence. 


Measured  Data  Model  1  - Model  3 


Frequency  (Hz) 


Figure  16.  Insertion  loss  of  0.635  cm  thick  PR- 1574  panel  with  low  (0.015)  loss 

tangent. 
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D.  NON-NORMAL  INCIDENCE  MODELS  AND  MEASUREMENTS 


Equation  13  (model  3)  is  capable  of  calculating  insertion  loss  for  non-normal 
incidence  cases  as  well.  As  in  section  C  above,  insertion  loss  measurements  were  made 
on  panels  of  several  material  types  at  prescribed  incidence  angles,  and  are  shown  below 
compared  to  experimental  measurements. 

The  case  of  the  Nylon  panel  at  non-normal  incidence  is  shown  in  Figure  17.  An 
insertion  loss  peak  at  160  kHz  in  Figure  17  is  significantly  overestimated  by  the  model, 
indicating  that  the  model  input  values  are  still  not  exact.  Since  the  normal  incidence 
model-to-measurement  match  was  good  and  hence  validated  the  values  used  for 
longitudinal  sound  speed  and  density,  this  spike  in  the  insertion  loss  data  likely  indicates 
the  Nylon  shear  estimates  require  update. 


Figure  17.  Insertion  loss  of  1.905cm  thick  Cast  Nylon  6  panel  at  15°  incidence. 


The  case  of  a  steel  plate  at  non-normal  incidence  is  shown  in  Figure  18.  Despite 
measurement  fluctuation  similar  to  that  seen  in  the  normal  incidence  case,  and  attributed 
to  the  same  cause,  the  match  between  the  model  and  the  measured  data  is  considered 
good,  indicating  shear  estimates  are  validated. 
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Figure  18.  Insertion  loss  of  1.0cm  thick  A36  carbon  steel  panel  at  20°  incidence. 

Once  again,  the  polyurethane  panel  measurement  results,  plotted  in  Figure  19, 
have  the  largest  deviation  from  model  3.  This  non-normal  incidence  case  shows  results 
similar  to  the  normal  incidence  case  and  further  supports  the  hypothesis  of  an  inaccurate 
attenuation  input. 
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Measured  Data - Model  1  - Model  3 
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Figure  19.  Insertion  loss  of  0.635cm  thick  PR1574  panel  at  30°  incidence. 
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IV. 


COMMENTS  AND  DISCUSSION 


A.  DESIGN  IMPROVEMENT 

Although  model  3  is  valid  for  multi-layered  materials,  no  actual  panels  of  such  are 
available  for  test  and  subsequent  model  comparison.  However,  a  comparison  of  several 
layered  models  provides  insight  into  the  effect  of  each  layer. 

A  benefit  of  using  a  layered  approach  is  the  ability  to  detennine  an  effective 
bending  stiffness.  Most  acoustic  windows  are  secured  to  the  hull  structure  such  that  the 
window  can  be  considered  a  plate  with  fixed  edges.  Structural  loads  carried  by  the 
window  limit  the  resulting  bending  displacement  criteria  to  ensure  geometric  stability  and 
equipment  clearances  are  maintained.  Displacement  and  bending  rigidity,  D,  are 
classically  related  by  the  plate  equation  given  in  Equation  24. 


p  =  V2D  V2w 
D  =  bending  rigidity 
w  =  displacement 
p  =  pressure  load 


Equation  24.  Classical  (Kirchoff)  plate  bending  equation 


N 

D  —  1/3  Z  (Ek  /  ( 1  -vk2))  (zk3  -  zk.,3) 

k=l 

k  =  Layer  number 

Equation  25.  Bending  rigidity  equation  for  plates 

Equation  25  is  based  on  the  work  by  Jones  [15]  which  uses  Kirchoff  s  thin  plate 
assumptions,  generally  considered  to  be  satisfied  for  plates  with  lateral  dimensions  at 
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least  20  times  the  thickness.4  The  bending  rigidity  of  a  plate  with  n  layers  can  be  broken 
into  n  separate  integrals  and  summed  to  determine  D  for  the  composite  section.  For 
example,  if  PR-1574  polyurethane,  1.905cm  thick,  is  sandwiched  between  two  .0635cm 
thick  Aluminum  plates,  the  effective  bending  rigidity  is  calculated  by  breaking  the 
integral  into  three  parts,  D1-D3,  and  summing  the  terms. 

For  this  case,  as  shown  in  Table  7  using  the  inputs  listed  in  Table  8,  D  for  the 
composite  section  increases  by  a  factor  of  -50  relative  to  a  single  layer  of  equal  thickness 
of  PR- 1574. 


D1 

4834.00 

D2 

151.63 

D3 

4834.00 

Djotal  for  Composite  Section 

9819.63 

D  Constant  section 

184.02 

Table  7.  Bending  rigidity  for  an  ALUM-PR1 574-ALUM  composite  section  and  a 

PR1574  constant  section  for  a  given  thickness. 


Material 

Density 

(g/cc) 

Young’s 

Modulus 

(GPa) 

Tan  8 

Layer 

Thickness 

(cm) 

Poisson ’s 
Ratio 

Layer  1 

2.7 

70 

.002 

.0635 

0.33 

Layer  2 

1.03 

0.20 

.1 

1.905 

0.49 

Layer  3 

2.7 

70 

.002 

.0635 

0.33 

Table  8.  Layered  material  1 : ALUM -PR1 574-ALUM. 


The  discussion  above  highlights  the  structural  benefits  to  a  layered  multi-media 
approach.  The  acoustic  perfonnance  can  be  evaluated  using  Table  8  inputs  for  a  3  layered 
system.  The  insertion  loss  is  calculated  for  various  angles  using  model  3  and  shown  in 
Figure  20. 


4  Plates  not  meeting  the  dimension  criteria  and/or  those  with  relatively  low  shear  rigidity  require 
‘thick’  plate  analysis  to  determine  bending  displacements.  See  Reissner  [16]. 
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Frequency  (Hz) 


Figure  20.  Model  3  insertion  loss  prediction  of  ALUM-PR1574-ALUM  ’sandwich’ 

composite  with  properties  as  shown  in  Table  8. 

Insertion  loss  of  the  ALUM-PR1574-ALUM  composite  section  shown  in  Figure 
20  increases  across  the  band  to  exceed  the  benchmark  of  2dB  at  only  a  small  portion  of 
the  high  end  of  the  band  for  the  angles  shown.  A  contour  plot  of  the  same  data,  shown  in 
Figure  21,  provides  a  better  display  of  the  data  and  reveals  that  this  benchmark  is  also 
violated  in  a  larger  part  of  the  frequency  band  at  ~  16°  angle  of  incidence.  The  sharp 
increase  in  insertion  loss  at  this  angle  can  be  shifted  in  angle  by  modifying  the  modulus 
of  layers  1  and  3,  as  evidenced  in  Figure  22  where  the  modulus  is  lowered  by  -20%. 
Lowering  layer  1  and  3  modulus  by  a  full  order  of  magnitude  (to  7GPa)  results  in  the 
complete  elimination  of  the  insertion  loss  peak,  as  seen  in  Figure  23.  This  suggests  that 
the  sharp  differences  in  modulus,  and  hence  acoustic  impedance,  between  the  outer  layers 
(1  and  3)  and  the  core  layer  (layer  2)  produce  internal  reflections  that  interfere  with  the 
sound  transmission  through  the  material. 
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Figure  2 1 .  Contour  plot  of  insertion  loss  (dB)  prediction  of  ALUM -PR1 574-ALUM 
'sandwich'  composite  with  properties  as  shown  in  Table  8. 
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Figure  22.  Contour  plot  of  insertion  loss  (dB)  prediction  of  ALUM-PR1 574-ALUM 
’sandwich’  composite  with  properties  as  shown  in  Table  8,  but  with  the  Young’s 
modulus  of  layers  1  &  3  changed  to  from  70  to  55GPa. 
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Figure  23.  Contour  plot  of  insertion  loss  (dB)  prediction  of  ALUM -PR1 574-ALUM 

'sandwich'  composite  with  properties  as  shown  in  Table  8,  but  with  the  Young’s 
modulus  of  layers  1  &  3  changed  to  from  70  to  7GPa. 

If  the  outer  layer  material  is  switched  from  Aluminum  to  A3  6  steel,  with 
properties  shown  in  Table  9,  D  increases  by  a  factor  of  -150,  relative  to  a  single 
equivalent  thickness  layer  of  PR- 1574.  Insertion  loss  values,  however,  increase  to 
unacceptable  levels,  as  shown  in  Figure  24. 
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Material 

Density 

(g/cc) 

Young’s 

Modulus 

(GPa) 

Tan  8 

Layer 

Thickness 

(cm) 

Poisson ’s 
Ratio 

Layer  1 

7.7 

207 

0 

.0635 

0.29 

Layer  2 

1.03 

0.20 

.1 

1.905 

0.49 

Layer  3 

7.7 

207 

0 

.0635 

0.29 

Table  9.  Layered  material  2:  Steel-PRl 574-Steel 
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Figure  24.  Contour  plot  of  insertion  loss  (dB)  prediction  of  STEEL-PR1 574-STEEL 
’sandwich’  composite  with  properties  as  shown  in  Table  9. 


These  two  examples  demonstrate  that  a  layered  system  can  be  designed  to  offer  a 
structurally  improved  product,  relative  to  the  mono-layer  design,  with  acceptable  acoustic 
performance.  In  addition  to  improvement  in  bending  rigidity,  designs  which  use  high 
strength  outer  layers  offer  improved  performance  relative  to  environmental  loads  as  well. 
For  example,  a  system  of  polyurethane  sandwiched  between  metallic  layers  (i.e.,  steel, 
Alum)  offers  greatly  reduced  water  absorption  by  the  polyurethane,  the  ability  to  accept 
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off-the-shelf  (OTS)  surface  treatments  for  fouling  prevention,  greater  resistance  to  in- 
service  wear/tear,  greatly  improved  chemical  resistance  and  reduced  design  restrictions  in 
mounting  schemes. 

B.  FUTURE  WORK 

The  next  logical  step  in  this  study  is  to  improve  the  quality  of  model  inputs  for  the 
layer  design  presented  in  Table  8  and  to  subsequently,  upon  acceptable  model  results, 
manufacture  a  layered  panel  according  to  that  design  and  test  for  both  acoustic  and 
structural  perfonnance.  This  step  results  in  model  verification  of  a  layered  system  after 
which  the  model  can  be  used  with  confidence  to  optimize  performance  for  application  in 
specific  cases. 

It  is  important  to  note  that  the  model  does  not  consider  the  compatibility  of  layers 
with  each  other.  Rather,  it  assumes  perfect  bonding  between  layers.  For  example,  the  case 
of  PR- 1574  and  Aluminum  is  known  to  be  a  pair  of  materials  that  can  be  bonded  together 
with  little  difficulty.  This  may  not  be  the  case  for  other  combinations  of  high  strength 
outer  layers  and  polyurethane  or  rubber  middle  layers.  Additionally,  there  likely  exists  a 
practical  limit  on  minimum  material  thickness. 
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APPENDIX.  MATLABCODE 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 


%%%%  This  program  models  the  transmission  of  sound  thru  a  seawater  bound  % 
%%%%  section  containing  layer(s)  of  isotropic  material  and  calculates  %%% 
%%%%  a  transmission  coefficient  and  an  insertion  loss  as  a  function  of  %% 
%%%%  frequency  and  angle  of  incidence  . 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 


clear 

%%%First  Input  is  the  number  of  layers 

numlayer=  input  ('  Enter  the  number  of  layers  in  the  solid  material:  '); 

numfreqsteps=  200;  %%%sets  the  resolution  of  the  frequency. ..this  is  the  number  of  rows  in  the 
IL  matrix 

numanglesteps=  200;%%%sets  the  resolution  of  the  angle  ...this  is  the  number  of  columns  in  the 
IL  matrix 

angleinc=.  00349; 

csw=1500;  %%%  Speed  of  sound  in  seawater 
cswT=0  ;  %%%  Shear  sound  speed  in  seawater 


%%%Units  used  throughout  are  Kilograms, Meters,  seconds, 


%%%%%%%%%  I N  PUTS 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%%%%%%%%% 

for  n=1:numlayer 


Ein(n)=input  (['  Enter  Tensile  storage  modulus  in  GPa  units  for  layer  #',num2str(n),':  ]); 

E(n)=Ein(n).*1e9; 

EtanD(n)=  input  (['  Enter  Tensile  tan  delta  for  layer  #',num2str(n),':  ]); 

Ei(n)=E(n)*(1+(j*EtanD(n))); 

v(n)=input  (['  Enter  Poisson  ratio  for  layer  #',num2str(n),':  ']); 


rhoi(n)=input  (['  Enter  Density  in  g/cc  for  layer  #  ',  num2str(n),  1); 
rho(n)=rhoi(n).*1000; 


thki=input  (['  Enter  thickness  in  cm  for  layer  #  ',  num2str(n),':  ]); 
thk(n)=  thki./l  00; 


end 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%%%%%%%%%%%%%%%%%% 


%%%%%%%%%%  START  CALCULATIONS  LOOP 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 


for  t=1  :numanglesteps 
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anglei(t)=angleinc.*t; 

deg(t)=anglei(t)*57.3;%%%1  rad  ~=57.3  degrees 
for  u=1  :numfreqsteps; 

frequency(u)=  50000+(u.*  200000./(numfreqsteps)); 
w=  2.*pi.*frequency(u); 
csw=1500; 
cswT=0  ; 

for  n=1:numlayer 


Pscale(n)=  v(n)./(  (1+v(n)).*(1-(2*v(n)))); 

lambdai(n)  =  Pscale(n).*Ei(n);  %%%%Lame's  constant 

lambda(n)=real(lambdai(n)); 

Pscale2(n)=(1./(1+v(n)))  *  ((v(n)/(1-(2.*v(n))))+  1); 
SpeedtermLi(n)=Pscale2(n).*Ei(n); 

%%%%%%Calculate  the  tan  d  for  the  compressional  wave  %%%%% 
LtanD(n)=  imag(SpeedtermLi(n))/real(SpeedtermLi(n)); 

%%  The  below  terms  are  used  in  several  locations  %%%% 

TD  =  1  +  ((EtanD(n))A2); 

TD1=sqrt(TD); 

%%%%%  Calculate  the  effective  Speedterm  (modulus)  for  the 
%%%%%  compressional  wave  using  the  undamped  modulus  and 
%%%%  the  Young'  loss  tangent%%%%%% 
SpeedtermL(n)=2.*real(Speedterml_i(n)).*(TD/(TD1+1)); 

SpeedL(n)=  sqrt((SpeedtermL(n))./rho(n)); 

XciL(n)=sqrt(  (rho(n)*(wA2)/(2*real(SpeedtermLi(n))))  *  ((TD1-1)/TD) ); 

%%%%%Calculate  the  effective  shear  wave  speed  using  the 
%%%%%Young's  modulus  and  poisson's  ratio 
Pscale3(n)=  1./(2.*(1+v(n))); 

Gi(n)=Pscale3(n).*Ei(n); 

SpeedtermT(n)  =(2*Gi(n))*  (TD/(TD1+  1)); 

SpeedT(n)=  sqrt(SpeedtermT(n)./rho(n)); 

XciG(n)=XciL(n);  %%  Shear  loss  =  Comp  loss  =  Young's  loss 
%%  due  to  simplification  made  in  Eq  17 


k(n)=  (w./SpeedL(n)); 
kappa(n)=  (w./SpeedT(n)); 

if  n>1  %%%%  If  there  is  more  than  one  layer  %%% 
angle(n)=asin  (SpeedL(n).*Snell(n-1 )); 
Snell(n)=sin(angle(n))./SpeedL(n); 
angleT(n)=  asin  (  SpeedT(n).*Snell(n)); 
else 

Snell  1=(sin(anglei(t)))./csw; 
angle(n)=asin  (SpeedL(n).*Snell1 ); 
Snell(n)=(sin(angle(n)))./SpeedL(n); 
angleT(n)=  asin  (  SpeedT(n).*Sneli(n)); 
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end 


M(n)=(k(n).*sin(angle(n)))-(j*XciL(n)); 

Pa(n)=(k(n).*cos(angle(n)))-(j*Xcil_(n)); 

Qa(n)=(kappa(n).*cos(angleT(n)))-(j*XciG(n)); 


Matrix  1  =  Potential  state  at  z=  thickness 

NOTE:  z  is  measured  in  the  direction  opposite  to  the  incident 
%%%%%%%%  propagation 
%%%%%%%%%%%%%%% 
%%%%%%%%%%%%%%%%%%%%%%%%%% 


z=thk(n); 

e1=exp((j*Pa(n)*z)-(XciL(n)*thk(n))); 

e2=exp(-j*Pa(n)*z); 

e3=exp((j*Qa(n)*z)-(XciL(n)*thk(n))); 

e4=exp(-j*Qa(n)*z); 


Kaya(n)=((((M(n)).A2)+((Pa(n)).A2)).*lambdai(n))+(2.*Gi(n).*((Pa(n))  A2)); 

A(1,1,n)=j*M(n).*e1; 

A(1,2,n)=j*M(n).*e2; 

A(1 ,3,n)=-j*Qa(n).*e3; 

A(1 ,4,n)=  j*Qa(n).*e4; 

A(2,1  ,n)=  j*Pa(n).*  el; 

A(2,2,n)=  -j*Pa(n).*  e2; 

A(2,3,n)=j*  M(n).*  e3; 

A(2,4,n)=  j*  M(n).*  e4; 

A(3,1  ,n)=-j*Kaya(n).*e1 ; 

A(3,2,n)=-j*Kaya(n).*e2; 

A(3,3,n)=-j*2.*Gi(n).*M(n).*Qa(n).*e3; 

A(3,4,n)=j*2.*Gi(n).*M(n).*Qa(n).*e4; 

Q2=(Qa(n)).A2; 

M2=(M(n))  A2; 

A(4,1  ,n)=-j*2*Gi(n).*M(n).*Pa(n).*e1 ; 

A(4,2,n)=  j*2*Gi(n).*M(n).*Pa(n).*e2; 

A(4,3,n)=-j*Gi(n).*(M2-Q2).*e3; 

A(4,4,n)=-j*Gi(n).*(M2-Q2).*e4; 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%%%%%%%%%%%%%%%%%%%%%%%%%% 

%%%%%%%%%%%  Matrix  2  =  Potentials  state  at  z=0  in  the  layer 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%%%%%%%%%%%%%%%%%%%%%%%%% 
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z=0; 

e1=exp((j*Pa(n)*z)-(XciL(n)*thk(n))); 

e2=exp(-j*Pa(n)*z); 

e3=exp((j*Qa(n)*z)-(XciL(n)*thk(n))); 

e4=exp(-j*Qa(n)*z); 


Kaya(n)=((((M(n)).A2)+((Pa(n)).A2)).*lambdai(n))+(2.*Gi(n).*((Pa(n)).A2)); 

A1(1,1,n)=j*M(n).*e1; 

A1(1,2,n)=j*M(n).*e2; 

A1  (1 ,3,n)=-j*Qa(n).*e3; 

A1(1,4,n)=  j*Qa(n).*e4; 

A1(2,1,n)=  j*Pa(n).*  el; 

A1(2,2,n)=-j*Pa(n).*e2; 

A1  (2,3,n)=  j*  M(n).*  e3; 

A1  (2,4,n)=  j*  M(n).*  e4; 

A1  (3,1  ,n)=-j*Kaya(n).*e1 ; 

A1(3,2,n)=-j*Kaya(n).*e2; 

A1(3,3,n)=-j*2*Gi(n).*M(n).*Qa(n).*e3; 

A1(3,4,n)=j*2.*Gi(n).*M(n).*Qa(n).*e4; 

Q2=(Qa(n)).A2; 

M2=(M(n)).A2; 

A1  (4,1  ,n)=-j*2*Gi(n).*M(n).*Pa(n).*e1 ; 

A1  (4,2,  n)=  j*2*Gi(n).*M(n).*Pa(n).*e2; 
A1(4,3,n)=-j*Gi(n).*(M2-Q2).*e3; 

A1(4,4,n)=-j*Gi(n).*(M2-Q2).*e4; 


B(:,:,n)  =A(:,:,n)*inv(A1(:,:,n)); 


if  n>1 

Bi(:,:,n)=B(:,:,n-1)*B(:,:,n); 

else 

Bi(:,:,n)=B(:,:,n); 

end 

end 

B1=Bi; 

B3=abs(B1); 

Y1  i=(((-1  ,/B1  (4,1  )).*B1  (4,2).*B1  (2,1  ))+B1  (2,2)); 

Y2i=(((-1  ./B1(4,1)).*B1  (4,3). *B  1(2,1  ))+B1  (2,3)); 

Y3i=(((-1  ./B1(4,1)).*B1  (4,2). *B  1(3,1  ))+B1  (3,2)); 

Y4i=(((-1  ./B1(4,1)).*B1  (4,3). *B  1(3,1  ))+B1  (3,3)); 

Pn=  (w/csw).*cos(anglei(t));%%%Pn  in  the  first  water  layer 
P1  =  (w/csw).*cos(asin  (csw.*Snell(numlayer)));%  PI  in  last  sw  layer 
Mn=  (w/csw).*sin(anglei(t));%%%Mn  in  the  first  water  layer 
rhon=1000; 
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rhol =1 000; 

lambdan  =rhon*(csw.A2); 

Kayn=((Mn.A2)+  (Pn.A2)).*lambdan; 

%%%Note:  Kayn  is  Kay  for  the  first  of  the  two  bounding  water 
%%%layers.  G  is  zero  in  water  so  the  expression  is  truncated 
%%%compared  to  Kaya(n)  used  above.  Lambda  is  found  from  the  sound 
%%%speed  in  water. 

H3=(-Kayn*((-Y1i.*P1)+(Y2i.*rho1*(wA2)) )); 

H4=(  Pn*((-Y3i.*P1  )+(Y4i.*rho1  *(wA2)))); 

Transl  =(-2.*rho1*Pn*Kayn); 

Trans2  =rhon.*(  (H3)+  (H4)); 


Trans  =  Trans1/Trans2; 
TranslL(u,t)=(abs(Trans)); 

IL  (u,t)=20*log10(1/TranslL(u,t)); 


end 

end 

%%%%%%%%%%%%%%%%%%%%%END  CALCULATIONS  LOOP 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 


%%%%%%%%%%%%%  PLOTS  %%%%%%%%%%%%%%%%%%%%%% 


subplot(1 ,3,1 ),  plot(frequency,IL) 
subplot(1 ,3,2),  imagesc(deg, frequency, IL) 
subplot(1 ,3,3),  surf(deg .frequency, IL) 
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